Phonon-assisted tunneling in an isolated double dot system 
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Phonon-assisted tunneling rates are evaluated for a well isolated double dot system defined in 
a GaAs semiconductor heterostructure of finite thickness. A separable model for the confining 
potential allows accurate determinations of doublet electron wavefunctions and energies. It is found 
that at small doublet energies the piezoelectric rates due to flexural modes give the dominant 
contribution. For small slab thicknesses the predicted rates are up to two orders of magnitude larger 
than for very thick slabs. 
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I. INTRODUCTION 

The problem we consider is the transfer of an electron 
between two dots separated by a tunneling barrier. The 
experimental situation we have is mind is that of two con- 
tiguous dots defined by electrostatic confinement, as in 
the experiments of references [1-5]. For open dots, where 
a current can be measured between leads connected to 
each dot, the problem has been well studied both experi- 
mentally and theoretically [6-8]. More recently however, 
experimental results have been obtained for an isolated 
double dot system in a GaAs heterostructure [3]. Trans- 
fer was forced by varying the gate voltages on one of 
the dots, and charge transfer detected by a non-invasive 
voltage probe. Since this system has obvious interest as 
a possible experimental realization of a q-bit and or as a 
half cell of a Quantum Cellular Automaton, a good un- 
derstanding of all mechanisms leading to charge transfer 
is desirable. 

The simplest mechanism for transfer is direct tunnel- 
ing, also called elastic tunneling: for a critical arrange- 
ment of the gate voltages some electron level has a wave- 
function equally spread over the two dots, and slight 
changes in the applied voltage suffice to confine the wave- 
function to one or the other of the two. For this to oc- 
cur it appears necessary to have a doublet of electron 
states degenerate in energy at the critical configuration. 
However, initial experimental studies already showed the 
importance of inelastic tunneling, mediated by acoustic 
phonons, when the energy separation between the dou- 
blet states is approaching the minimum [1]. Although 
the present generation of experiments has been carried 
out on two dimensional electron gases (2DEG) defined 
in very thick heterostructures, similar experiments may 
be done in future [9] on 2DEG defined in heterostruc- 
tures consisting of much thinner slabs. It is of interest 
to examine the effect that changes in the spectrum of 
phonon modes, due to finite slab thickness, would have 
on the phonon-assisted tunneling rates. A study along 
these lines has recently appeared for open double dots 
[7]. The purpose of the present work is complementary: 



we will address the problem of a strongly isolated double 
dot system using a more realistic model for the confining 
potential. 

II. TWO COUPLED DOTS: SEPARABLE MODEL 

In an earlier study [10] on a device similar to that in 
the experiment of Gardelis et al [3], we found that to 
good approximation a strongly isolated dot was well de- 
scribed by a confining potential of roughly rectangular 
shape and constant depth surrounded by very steep bar- 
riers. To account for leakage of electrons through one 
side of the dot it was sufficient to use a barrier of finite 
height and width derived from Poisson Thomas-Fermi 
simulations. Guided by the success of that model in ex- 
plaining the lifetimes of the emitted electrons, we will use 
here a similar confining potential. For simplicity and to 
achieve satisfactory numerical accuracy, we have chosen 
a rectangular shape for the barrier separating the two 
dots. Therefore, in the plane of the 2DEG we consider 
a double dot of rectangular shape, and define our x and 
y axes along the directions of its sides. We further set 
the confining potential, V(x, y), to be infinite outside the 
boundaries of the rectangle and inside to depend only on 
x: the profile V(x) is shown in Fig. 1. Dot A on the left 
has sides wi x l y , and dot B on the right, w r x l y . The 
square barrier separating them runs from —Wb/2 to u>&/2. 
For convenience we set the zero of energy at the top of 
the barrier, so that inside dots A and B the potential 
is Vd < 0. The electron wavefunctions, in the envelope 
function approximation, ^>{x,y,z) — (j)(x)ip(y)x{z), are 
obtained by solving a 3D separable Schrodinger equation. 

For simplicity, we chose the well known Fang-Howard 
[12] parametrization for x( z )i with parameters appropri- 
ate for the Cavendish experiment. We give its explicit 
form in Section III. Due to the simple form of the con- 
fining potential along y, the tp(y) are the well known 
eigenf unctions of the infinite square well. Along the x 
axis: 
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b"{x) + V{x)<P{x) =E^<f>(x) 



(1) 



and V(x) = Vd when — wi — Wb/2 < x < —Wb/2 
and Wft/2 < x < Wb/2 + w r , whereas V(x) = when 
—wt,/2 < x < Wb/2. The obvious analytic forms of 
4>{x) in the "zones" (A, B or barrier) arc matched at 
the boundaries (wavefunction and derivative) leading to 
explicit equations for the energy that are easily solved 
numerically to the required accuracy, even for sizeable 
barriers like those studied here. Once these eigenener- 
gies have been determined it is a simple matter to find 
the fraction of the normalization of the wavefunctions in 
each dot. The assignment of each solution to a level in 
dot A or B is in most cases unambiguously determined by 
the fractions of normalization. By counting the number 
of nodes in the corresponding dot we assign to each solu- 
tion an effective quantum number n; (left) or n r (right.) 
The level crossings which occur for special values of the 
parameters are detected by the jumps, from almost zero 
to almost one, of those normalization fractions. There is 
always negligible normalization inside the barrier. 
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FIG. 1. Confining potential profile along the x direction. 

The energy of a given level is the sum of the eigenval- 
ues corresponding to the three eigenfunctions into which 
^(x, y, z) is factorized. However, in a doublet, the differ- 
ence of energy between initial and final states has equal 
contributions from the y and z components since they 
must have the same tp(y) and x( z )- Therefore in what 
follows E if = E { ;l a} -E%l al . 

For definitencss, in the present study we will use the 
width of dot B, w r , as the parameter whose variation 
simulates the changes induced by the variation of one of 
the confining gate voltages. In Fig. 2 we show the level 
crossing chosen to present the results of our studies. The 
parameters for the double dot configuration are: Vd = 
—35 meV, wi = 404 nm, l y — 350 nm, Wb — 24 nm and 
w r is varied around 380 nm. These parameter values are 
similar to those used in previous simulations of double 



dot devices also measured at the Cavendish [10,11]. The 
doublet shown corresponds to n; = 17, n r — 16 and 
with tp(y) = \j2jl y cos(ny /l y ), has a total energy close 
to the estimated Fermi level in the double dot system. 
The continuous lines in the figure are the doublet levels, 
whereas the dashed lines are inserted to guide the eye 
and would correspond to the "unperturbed values" of an 
ideal uncoupled double dot system. 

In a simple two-level system described by a hamilto- 
nian: 
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AE = \/(E A — E B ) 2 + 4t 2 is the separation of the two 
eigenstates. When the two unperturbed values cross, 
Ea = Eb, and t — AE/2. This allows one to determine 
an effective tunneling width which characterizes the cou- 
pling between the two dots. From the data shown in Fig. 
2 we find t = 2.2 ^eV. 
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FIG. 2. Energies. Doublet avoided level crossing, n r 
ni = 17. 
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To simulate an experiment we will assume for instance 
that the initial conditions correspond to w r — 380.4 nm 
and that the lowest level in the right part of figure 2 (cor- 
responding to dot B) is occupied by one electron and the 
upper level in the same right part of the figure (dot A) is 
empty. We then reduce the value of w r , shrinking the size 
of dot B but leaving dot A unchanged. The energy of the 
lowest level increases linearly as w r decreases so long as 
the occupation of the corresponding wavefunction stays 
close to unity in dot A, and stabilizes when the electron 
tunnels to dot B at about w r — 380 nm. This is direct 
or elastic tunneling. Inelastic tunneling due to phonon 
absorption will take place when w r is still above 380 nm, 
via transitions from the lower level (electron in dot B) 
to the upper level (electron in dot A.) In our model we 
are not dealing explicitly with the coulomb energies of 
the electron in each well. For the configurations where 
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the rates are maximal the initial and final electron wave- 
functions are equally spread over the two dots, so we can 
assume that the contribution of Coulomb energy to the 
doublet separation cancels. 



III. ACOUSTIC PHONON RATES FOR A FINITE 
SLAB HETEROSTRUCTURE 

According to Fermi's golden rule, the rate for phonon 
assisted tunneling is: 

w = — | < $ F * F \H e _ ph \$* > | 2 S(E -E f t M (3) 

where the bras and kets indicate initial and final states of 
the electron (\^f >) plus phonons (|$ >) and the =F sign 
applies to phonon emission or absorption. Since GaAs is 
piezoelectric we will consider both the deformation po- 
tential and the piezoelectric contributions to the electron- 
phonon Hamiltonian. The specific forms of these rates for 
a finite slab have been thoroughly discussed in the work 
of Bannov et al. [13] which we follow closely. The hamil- 
tonians describing both contributions can be written in 
second quantized form: 
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with the superscripts d and p indicating deformation and 
piezoelectric contributions and the subscripts d, f and 
s dilatational, flexural and shear modes. The z-axis is 
perpendicular to the plane of the 2DEG and the faces 
of the slab, and the origin is in the middle of the slab. 
Vectors in the x, y plane are indicated by the subindex ||, 
and n is the mode index. The expressions for the various 
T's are found in equations 33-34 (deformation potential) 
and 39-41 (piezoelectric) of [13]. Particularizing to the 
separable electron wavefunctions of our model: 
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J2 -Mil \<x(z)\T%f /s (q=,n,z)\ X (z)> 



x (N(hw n ) + l±l) S{E-E't Kuj n ) , (5) 

where we have evaluated the matrix elements correspond- 
ing to the phonon occupations, and defined 
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<M^M2/)|e^- F 1<M*Wy) > 



(6) 



We remark that Bannov et al. [13] chose the axes in the 
2DEG plane so that q\\ ~ (q x , 0), whereas we have aligned 
the x and y directions along the sides of the double dot. 
Therefore, the q x appearing in their work correspond to 
|aj|| = </|| in ours. Accordingly, we write: 



M\\ = | < (j) F (x)\e^ xcos6 ^(j}{x) > | 2 

x | < i>(y)\e lq u vsin6 m>(y) > | 2 



(7) 



To show the coupling constants explicitly, we rewrite the 
T's as follows: 
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with explicit expressions for the 7's to be found in the 
above mentioned equations of [13]. The normalization 
factors, F d ff/ S are given in their appendix. In this revised 
notation, the deformation rate becomes 



x \<x{z)\l^\z)\ X {z)>\ 2 



x (N(hu n ) + \ t ^ S(E-E fT hw n ) , (9) 

where we have omitted the subindexes d/f that appear 
each time a mode n is referenced. A similar expression 
applies to the piezoelectric rates. Furthermore, for an 
infinite slab area, one has to replace 1 /A by the in- 
tegral (1/(2tt) 2 / dq\\ = (1/(2tt) 2 / q\\ dq\\ d6\\, so that 
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We then arrive at 
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To determine the expectation values of the 7's we have 
chosen a Fang-Howard wavefunction to represent x( z )' 



= , z < z . 



z > z 



(13) 

with z = a/2 — z and z is the depth of the 2DEG referred 
to the surface of the slab. We have chosen z = 70 nm 
and 6=1/4 nm -1 . 

We solved the algebraic equations (13), (14) and (18) of 
[13] numerically, so as to obtain the dispersion relation 
?iu) n {q\\) and corresponding expressions for the other mo- 
menta q t and qi appearing in the explicit expressions for 
the ^' p \z) and the normalization factors F 2 . 
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IV. RESULTS 



preserved by less schematic double dot potential profiles. 



Before presenting results for the rates, it is interesting 
to discuss some general features of the terms appearing in 
cq. 1 1 and their dependence on the doublet energy split- 
ting, Eif. For a given mode, the corresponding disper- 
sion relation Eif = tuv n (q\\;a) determines the momenta 
q\\, qt and qi. As shown in figure 2, to a good approxi- 
mation one finds the same value of the splitting energy 
for pairs of points located symmetrically with respect to 
the crossing point at w r ^ cross ~ 380 nm. One such pair 
is w r ,i — 379.8 nm , w r .2 = 380.2 nm. For these two the 
values of gy are nearly the same. This implies that all fac- 
tors contributing to eq. 11 except My will also be almost 
equal since they depend only on gy . The matrix element 
.My, involves the doublet wavefunctions <f>(x) and (/>f(x) 
and these depend on the value of w r . It turns out how- 
ever that due to the simplicity of the potential profile 
chosen for V(x), the products of wavefunctions for the 
paired doublets are also very similar. This is shown in 
Fig. 3, where we have plotted the overlaps <f>F{x)<j){x) 
for the doublets w r _i — 379.8 nm ( continuous line ) and 
w r ,2 = 380.2 nm (dashed line). For clarity of the figure 
we only show the range —100 < x < 100 nm, but the 
two lines practically coincide in most of the range of x's. 
If we had chosen uv,i = 379.77 nm, so as to make the 
two Eif even closer, the two lines in the plot would be 
indistinguishable 

0.0015 | , , , 1 
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FIG. 3. Overlaps of wavefunctions. Continuous line: the 
<Pf{x)4>(x) of the pair corresponding to w r = 379.8 nm. 
Dashed line: same for the pair at w r — 380.2 nm. 

The result of this symmetry is that the matrix clement 
.My, is nearly identical for the two cases, since it is deter- 
mined by the overlaps. In summary: the rates are very 
similar for pairs with w r symmetric with respect to the 
crossing point. Part of this is due to the simplicity of the 
potential profile chosen, but as long as the potential pro- 
file at the barrier is symmetric and each dot has a wide 
flat potential well, we expect that this symmetry will be 
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FIG. 4. Matrix element A4u for various barrier widths: 
Continuous line w b — 24 nm (t e // = 2.2/xeV) , dotted line: 
w b = 18 nm (t e ff = 7.6/ieV), dashed line: w b — 36 nm 
(t e/ / = 0.2/xeV). 

Fig. 4 shows .My for three barrier widths. One sees 
that indeed the curves are symmetric with respect to 
Eif = 0, and that the magnitude depends strongly on 
Wb- However, their shape depends little on the barrier 
width, showing the expected increase at small doublet 
energies. To analyze the latter, it is useful to consider 
the contributions to .My for very low q\\. When gy — > 0, 
eq. 7 can be approximated by 

A4y ~ g 2 cos 2 0„ | < <j> F {x)\x\<j>{x) > | 2 . (14) 

Fig. 5 shows these matrix elements M. for a configu- 
ration with w r = 379.9 nm and tp(y) the lowest energy 
solution for the infinite square well. It can be seen that 
gy ~ 4.5 . 10~ 3 nm -1 (corresponding to Eif ~ 6.1/xeV) 
is not sufficiently close to and therefore eq. 14 gives 
only a rough approximation for most values of fly. For 
other values of w r the agreement is similar. Eq. 14 is 
nevertheless useful because it correlates M. y with the ex- 
pectation value of x, and the behaviour of the latter is 
easier to visualize: 

For large doublet splitting the overlap of the wavefunc- 
tions <fi(x) and <j)F{x) is small: one is mainly located in 
dot A, the other in B. As Eif decreases each wavefunc- 
tion begins to spread into the other dot. At the mini- 
mal doublet splitting, and assuming for sake of argument 
equal dot sizes, the wavefunctions become a symmetric 
and antisymmetric pair each with equal occupations in 
the two dots. The overlap is then maximal, but the prod- 
uct 4>{x)4>f(x) is of positive sign in one dot and negative 
in the other so the resulting overlap integral vanishes. In- 
troducing an x in the integration, as in eq. 14, makes the 
contributions from the two dots equal and therefore max- 
imizes the matrix element when the two eigenfunctions 
are maximally spread over the two dots. 
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FIG. 5. The matrix elements entering M v.s. Q\\ ; two dot 
system with w r = 379.9 nm. Continuous line: Full M. Close 
by dashed line: | < <f> F \e iq W cos0 \\ x \4> > | 2 . Short dashed line 
at top: | < Tp F \e tq \\ sine \\ y \ip > \ 2 . Dotted line: approximate 
M of eq. 14. 
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FIG. 6. Matrix element Af as a function of qu for the lowest 
flexural mode. Continuous line: exact numerical calculation. 
Dashed line: small q\\ expansion. 

To complete the analysis, we show in Fig. 6 the matrix 

element : Af = \Ff < x{ z )\ln\z)\x{z) > |, for the lowest 
flexural mode. As will be shown in next section this 
mode gives the dominant contribution to the rate at small 
doublet splittings. At very small q\\ one can perform 
power series expansions of the quantities entering into 
this matrix element and, neglecting the 2DEG thickness, 
one finds 

\r jjff cosh (gn Q / 2 ) cosh (g||( a / 2 - go)) n ^ 
^ (i-( St / S; ) 2 ) 2 ^ ' [ ' 

where si and s t are the velocities of bulk longitudinal 
and transverse waves. This approximation is shown as a 



dashed line in Fig. 6 and reveals the cause of the sharp 
increase when q|| — > 0. In conclusion, the two matrix 
elements, M\\ and Af contribute to the strong increase 
of the piezoelectric rates at small doublet energies, to be 
shown next. 

Rates: Following ref. [3], we have chosen T = 100 mK 
for all our rate calculations. Also, for defmiteness, we 
have assumed that one of the symmetry axes of GaAs is 
parallel to z, whereas the other two are in the plane of the 
2DEG, rotated by an angle 8 sym = ir/5 with respect to 
the x, y axes introduced above. The material parameters 
and coupling constants for GaAs were taken from [14,15]. 
As discussed at the beginning of this section, to the accu- 
racy of the logarithmic graphs to be presented, the rates 
for pairs of states with the same energy difference -E 1 , / are 
indistinguishable. To present emission and absorption 
rates in the same graph, we adopt the convention that 
negative (positive) Ei f correspond to phonon absorption 
(emission). Again, to the accuracy of the graphs, the 
difference between the two rates for the same \Eif\ is 
entirely due to the factor N(E if ) + 1/2 ± 1/2. 

In Figs. 7 to 9 we present the rates for a slab of thick- 
ness a = 500 nm, for the doublet shown in Fig. 2. The 
2DEG is located at a depth zq = 70 nm below the surface 
and ip(y) is the lowest eigenfunction of the square well. 
Note the different scales in Fig. 7 from the rest. This 
reflects the dominant role of piezolelectric rates due to 
flexural modes, particularly at small doublet energies. In 
Fig. 7, we show the piezoelectric rates due to the flexu- 
ral modes. The various dashed and dotted lines show the 
separate contributions for each mode, and the continuous 
line is the sum of all. Note that only the lowest flexural 
mode has a dispersion relation such that hui — > when 
<7|| — > so that it contributes to all doublet separations. 
Higher modes tend to a finite energy when q\\ — > and 
therefore their contributions start at increasing non-zero 
threshold values of Eif. Note that at small \Eif\ the long 
dashed line corresponding to the lowest flexural mode is 
indistinguishable from the sum due to the negligible con- 
tribution of the remaining modes. 

At the beginning of this section we have already ex- 
plained the strong increase in this rate when Eif — ► 0. As 
the doublet separation increases the total rate decreases, 
with oscillations that only appear smooth due to use of a 
logarithmic plot, but which have significant features: at 
around 30 /ieV one sees a sharp jump due to the threshold 
of the third flexural mode and a pronounced dip is seen 
at 60 /LteV due to the oscillating contribution of the fifth 
mode. For the same flexural modes, the rates at small 
doublet energies are roughly eight orders of magnitude 
larger than those of the deformation potential interac- 
tion. Only at the threshold ~ 0.1 meV does the defor- 
mation rate reach a value comparable to the piezoelectric 
rate. The deformation rates might perhaps be more rele- 
vant in dots defined on a less piezoelectric semiconductor, 
such as In As. 

Figs. 8 and 9 show the piezoelectric rates due to shear 
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and dilatational modes. They are of similar magnitude 
to each other, and only beyond 50^teV do they become 
comparable to the flexural ones. Note also that they be- 
come very small at the lowest doublet separations. The 
deformation rates due to the dilatational modes are al- 
ways one to two orders of magnitude smaller than the 
corresponding piezoelectric rates. 

The dependence of the rates on slab thickness is shown 
in Figs. 10 to 12. In Fig. 10 we show piezoelectric rates 
due to flexural modes, for a fixed doublet splitting energy 
that corresponds to w r = 379.9 nm (Eij ~ 6.1 /ueV.) The 
contribution from the lowest flexural mode (continous 
line) dominates the others (dashed and dotted lines) over 
the whole range of thicknesses considered. Interestingly, 
increasing a from 500 to 2000 nm reduces the total rate 
by two orders of magnitude, showing the importance of 
finite slab thickness in determining rates. 

Fig. 11 examines the dependence of that reduction on 
the 2DEG depth. Bringing zo closer to the surface in- 
creases the rates further, by about 25% in going from a 
depth of 70 nm to 35 nm . Finally, Fig. 12 shows the dif- 
ference in rates between a "surface" and a "bulk" 2DEG: 
the location of the 2DEG has been set at zq — a/2 — 180 
nm and a has been varied. The configuration, a = 500 
nm, which sets zq = 70 nm, was our standard configura- 
tion in all previous graphs. Now, as a is increased, the 
2DEG stays near the center of the slab, and its distance 
to the surface increases. One can see in Fig. 12 two 
distinct regimes: "surface": for small slab thicknesses 
the total rate due to all flexural modes decreases expo- 
nentially by as much as five orders of magnitude ; ii) 
"bulk" beyond a ~ 2000 nm the rate stabilizes in or- 
der of magnitude, but oscillations of up to a factor of 2 
still persist. The contributions to the rate due to shear 
and dilatational modes (not shown) are totally negligible 
at small a and show a similar stabilization for a above 
2000 nm, where they are still one order of magnitude 
smaller than the flexural ones. Although experiments 
with a 2DEG at such depths are unrealistic, these results 
highlight the danger of estimating the rates for a double 
dot system without considering the influence of surface 
boundary conditions on the phonon spectrum. 

V. CONCLUSIONS 

We have studied the main features of the tunneling 
rates between a doublet of states in a well-isolated two 
dot system defined on a slab of finite thickness. At low 
temperatures, T = 100 mK, and for a piezoelectric ma- 
terial (GaAs), our results confirm the dominance of the 
piezoelectric over the deformation rates. At small dou- 
blet separations we find that the lowest piezoelectric flex- 
ural mode gives the dominant contribution. Keeping the 
2DEG at a fixed depth from the surface and increasing 
the slab thickness we find that the rates vary as much as 
two orders of magnitude, increasing as the slab thickness 



is decreased. Therefore in thin slabs the effective electron 
to acoustic phonon coupling is greatly enhanced. Com- 
plementarily, as a way to stress the relevance of bound- 
ary conditions on the phonon spectrum, when we keep 
the 2DEG near the middle of the slab and increase its 
thickness, we find that the "bulk" configuration leads to 
substantially smaller rates than the "surface" one. These 
findings should be useful in the design of future experi- 
ments on double dot systems supported on finite slabs. 
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FIG. 7. Piezoelectric rates due to flexural modes, as a func- 
tion of doublet splitting. Continuous line sum of all contri- 
butions. Dashed lines : separate contributions of each of the 
modes 



FIG. 10. Piezoelectric rates due to flexural modes, for fixed 
w r — 379.9 nm ( Eif = 6.1/ieV ) , as a function of slab 
thickness a. Continuous line: contribution of the lowest mode. 
Dashed lines : other modes. 
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FIG. 8. Piezoelectric rates due to shear modes, as a func- 
tion of doublet splitting. 
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FIG. 11. Piezoelectric rates for the lowest flexural mode, 
for fixed w r = 379.9 nm (Eif = 6.1/ieV), as a function of 
slab thickness a; z = 70 nm (continuous line), zo = 35 nm 
(dashed) . 
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FIG. 9. Piezoelectric rates due to dilatational modes, as a 
function of doublet splitting. 
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FIG. 12. Piezoelectric rates due to the flexural modes, for 
w r — 379.9 nm (Eif — 6.1/ieV) and zo = a/2 — 180 nm, 
against a. Dashed lines: individual modes; continuous line, 
total. 
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